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Abstract 

We study an asymptotic behavior of a special correlator known as the Emptiness Formation Probability 
(EFP) for the one-dimensional anisotropic XY spin-1/2 chain in a transverse magnetic field. This correlator 
is essentially the probability of formation of a ferromagnetic string of length n in the antiferromagnctic 
ground state of the chain and plays an important role in the theory of integrable models. For the XY 
Spin Chain, the correlator can be expressed as the determinant of a Toeplitz matrix and its asymptotical 
behaviors for n — ► oo throughout the phase diagram are obtained using known theorems and conjectures on 
Toeplitz determinants. We find that the decay is exponential everywhere in the phase diagram of the XY 
model except on the critical lines, i.e. where the spectrum is gapless. In these cases, a power-law prefactor 
with a universal exponent arises in addition to an exponential or Gaussian decay. The latter Gaussian 
behavior holds on the critical line corresponding to the isotropic XY model, while at the critical value of the 
magnetic field the EFP decays exponentially. At small anisotropy one has a crossover from the Gaussian to 
the exponential behavior. We study this crossover using the bosonization approach. 
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I. INTRODUCTION 

Although the Bethe Ansatz gives us important information about the ground state 

and thermodynamics of quantum integrable one-dimensional models, the calculation of correlation 
functions for these models is still an open problem. Considerable progress, however, has been made 
recently in this field (see P, 0] and references therein) . 

Two main approaches express correlation functions in integrable models in terms of determinants 
of Fr e d a„ lm operators flff or. ™, ipl e ^ Q. T.ese expressions are exact but are v ery 
complex. 

It was argued 0] that the simplest of the correlators in some integrable models is the so-called 
"Emptiness Formation Probability" (EFP). For an XYZ spin chain in a magnetic field which is 
defined as 

N 

i=l i 

the Emptiness Formation Probability represents the "Probability of Formation of Ferromagnetic 
Strings" 

P(n) = il> je-f n^|> ( 2 ) 
where Z = Tr {e - ^ 7 "}. At zero temperature T = it becomes 

P(n) = (0|n^H0>, (3) 

i=l 

that is, the probability that n consecutive spin sites are all found aligned downward in the ground 
state |0). 

It is conceivable that the study of this simple correlator will bring insights helpful to carry on 
the investigation of other correlators. But despite the fact that the EFP is the simplest of the 
correlators and despite considerable efforts devoted over the years to its study, there is still no easy 

3 



recipe for its calculation (see, for instance, 



12J). However, some progress has been achieved 



in finding an asymptotic behavior of P(n) at large distances n — * 00. 

For the XXZ spin chain in zero magnetic field (J x = J y , h = in ([T])). the EFP is found to 
be Gaussian P(n) ~ e~ an as n — ► 00 in the critical regime | J 2 I I J x I ctt zero temperature and 



exponential e ^ n at finite temperature ( 



1$, 



A qualitative argument in favor of Gaussian decay was given in Ref. |f| within a field theory 
approach. It was argued there that the asymptotics of the EFP are defined by the action of an 



optimal fluctuation (instanton) corresponding to the EFP. In the critical model, this fluctuation will 
have a form of a "n x n" droplet in space-time with the area A ~ n 2 and the corresponding action 
S ~ an 2 which gives the decay P(n) ~ e~ an2 . Similarly, at finite temperature the droplet becomes 
rectangular (one dimension n is replaced by an inverse temperature T -1 ) and the action cost is 
proportional to n, giving P(n) ~ e~^ n . This argument is based on the criticality of the theory 1 
and it is interesting to consider whether it could be extended to a non-critical theory. A naive 
extension of the argument would give the optimal fluctuation with space-time dimensions n x £ 
where £ is a typical correlation length (in time) of the theory. This would result in P(n) ~ e~^ n 
for non-critical theories, similarly to the case of finite temperature in critical regime. The rate of 
decay f3 would be proportional to the correlation length of the theory. 

In this paper we examine the relation between the asymptotic behavior of the EFP and criticality 
using the example of the Anisotropic XY spin- 1/2 chain in a transverse magnetic field 



N 



i=l 



H = E ( ^ k + ( V ) * -*E* W 



2 



N 



i=l 



where erf, with a = x,y, z, are the Pauli matrices which describe spin operators on the i-th lattice 
site of the spin chain and, for definiteness, we require periodic boundary conditions: of = &f +N 
(N >> 1). This model is, probably, the simplest model that has both critical and non-critical 
regimes in its 7 — h phase diagram. 

Using the mapping of the model (J3J to free one-dimensional fermions (Jordan- Wigner transfor- 
mation) one can express the EFP P{n) of this system in terms of a determinant of a n x n Toeplitz 
matrix. The asymptotic behavior of these matrices can be found using known theorems from the 

n 

theory of Toeplitz determinants. This approach and technique is exactly the one of Ref. |13|], 

where it was used to find the spin-spin correlation functions for the model @. This technique was 

n 

also used in Ref. [3] for the EFP in the case of the Isotropic XY model (Eq. @ with 7 = 0). In 
the latter work it was shown that the EFP decays in a Gaussian way for the critical theory (7 = 0, 
— 1 < h < 1). This case corresponds to one of the two critical lines in the 7 — h phase diagram 
of the model The other line is the critical magnetization line(s) (h = ±1). In the rest of the 
two-dimensional 7 — h phase diagram, the model is non-critical. 

We obtain that the EFP is asymptotically exponential in most of the phase diagram according 
to the naive expectations and that it is Gaussian only at 7 = in agreement with previous 
works on XXZ spin chains and Ref. 7]. However, on the critical lines h = ±1, in addition to 



1 More precisely, on the assumption that temporal and spatial dimensions of an instanton scale similarly. 
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the exponential decay, a pre-exponential power-law factor arises, with a universal exponent. The 
power-law prefactor is present in the isotropic case as well, but with a different exponent. 

The paper is organized in the following way: we briefly review the phase diagram of the XY 
Spin-1/2 Chain in Section [H] In Section IlIII we explain how one can express the EFP as the 
determinant of a Toeplitz matrix and review our results so that readers who are not interested 
in derivations can skip the next sections. In Section HVI we analyze the exponential decay of the 
EFP for the non-critical and critical phases of the anisotropic XY Model. In Section we derive 
in detail the asymptotic behaviors, including the pre-exponential factors, of both non-critical and 
critical parts of the phase diagram. In Section IVII we study a special line of the phase diagram 
on which the ground state is known exactly and compare the explicit results one can obtain using 
the exact ground state with the asymptotes of the EFP we derived in the previous sections. In 
Section IV11I we report on the already known results for the EFP of the isotropic XY model |2j • In 
Section [Villi we make contact with Ref. Q] using a bosonization approach to discuss the crossover 
as a function of n from the Gaussian to the exponential behavior of EFP for the case of small 
anisotropy 7. Finally, Section IIXI will summarize our results. For the reader's convenience we 
collect some results on asymptotic behavior of Toeplitz determinants which are extensively used in 
the rest of the paper in the first appendix. The second appendix presents the analysis of the finite 
temperature behavior of the EFP, which gives an expected exponential decay. The third appendix 
gives some mathematical details on the calculation of the stationary action in the bosonization 
approach of Section [Villi 

I I 

Some of the results presented in this paper have been announced in a previous letter without 
details, which will be given here. 



II. THE ANISOTROPIC XY MODEL 

The XY spin model defined by 0) has been solved in Il5|] in the case of zero magnetic field and 

Q n 
in the presence of a magnetic field. We follow the standard prescription 15] and reformulate 

the Hamiltonian Q) in terms of spinless fermions ifii by means of a Jordan- Wigner transformation: 



a+ = rp^k^k = ^ "Q f 2 ^ k _ A ; (5) 



k<j 



i 



2^ f Vi - 1, (6) 



3 



where, as usual, = (a x ± ia y )/2: 

N 

H = J2 (4^+i + + 7 44+i + 7 - 2/1 4^) ■ ( 7 ) 

i=l 

In Fourier components ipj = ^ ipq^ > we have: 

= ^ [2 (cosg - h) ip\ip q + ry sing ip\4- q ~ h sing ip- q ip q ■ (8) 
The Bogoliubov transformation 



■d 

X 9 = cos-^ ip q + i sin^ ip1 q (9) 



with "rotation angle" $ q 



e t ' dq = —(cos q - h + Z7 sing), (10) 



brings the Hamiltonian (JHJ) to the diagonal form ^gXgXg with the quasiparticle spectrum 



e q = y (cos g — /i) 2 + 7 2 sin 2 g. (11) 

We recognize from (|11|) that the theory is critical, i.e. gapless, for h = ±1 or for 7 = and \h\ < 1. 

In Fig. n w e show the phase diagram of the XY model, which has obvious symmetries 7 — > —7 
and h — > —h. However, the latter is not a symmetry of the EFP. Therefore, we show only the 
part of the diagram corresponding to 7 > 0. The phase diagram has both critical and non-critical 
regimes. Three critical lines fio (Isotropic XY model: 7 = 0, \h\ < 1) and Q± (critical magnetic 
field: h = ±1) divide the phase diagram into three non-critical domains, S_, So, and (h < — 1, 
— 1 < h < 1, and h > 1 respectively). Fig. ^also shows the line 7 = 1 (r^) corresponding to the 
Ising model in transverse magnetic field and the line 7 2 + h? = 1 (r p) o n which the wave function 
of the ground state is factorized into a product of single spin states lj| (we will consider this line 
in detail in Section IVTI) . 



The fermionic correlators are easy to obtain from (|8J). In the thermodynamic limit they read 
F jk = Wrfk) = -i(44 k ) = 1^ \ sin (12) 



4 



h 



i 



£2 



-i 



i 






\ 

\ 

\ 

\ 

\ 

\ 




1 


1 







FIG. 1: Phase diagram of the XY Model (only the part 7 > is shown). The theory is critical for h = ±1 
(f2±) and for 7 = and \h\ < 1 (Slo) • The line Ti represents the Ising Model in transverse field. On the line 
the ground state of the theory is a product of single spin states. 
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III. EMPTINESS FORMATION PROBABILITY IN THE XY MODEL 



We are mostly interested in the value of the correlator © at zero temperature (the non-zero 
temperature case is deferred to the second appendix). In terms of spinless fermions, one can express 
the EFP (j3J) as the expectation value over the ground state of the theory [?J 

n 

?W = (II# ( 14 ) 

3=1 

This expression projects the ground state on a configuration without particles on a string of length 
n and gives the meaning to the name "Emptiness Formation Probability". 

Let us introduce the 2n x 2n skew-symmetric matrix M of correlation functions 

/-IF G\ 

M = , (15) 

V-g if; 

where F and G are n x n matrices with matrix elements given by Fjk and Gjk from H12I13[) 
respectively. Then, using Wick's theorem on the r.h.s of (|14j) . we obtain 



P(n) = Pf(M). (16) 



n 

The Pfaffian 16] is defined as 



Pf(M) = ^(-1) P M P1P2 M P3P4 . . . M p2n _ lP2n , (17) 
p 

where P = {pi,p%, ■ ■ ■ ,P2n} is a permutation of {1, 2, . . . , 2n}, the sum is performed over all 
possible permutations, and (— 1)^ is the parity of the permutation. Using one of the properties of 
the Pfaffian we have 

P{n) = Pf (M) = v/det(M). (18) 
We perform a unitary transformation 

, f S A if 1 -A 

M' = UMUt = , U = -= , (19) 

where I is a unit n x n matrix and S n = G + iF and S n ^ = G — iF. This allows us to calculate 
the determinant of M as 

det(M) = det(M') = det(S n ) ■ det(S n f ) = |det(S n )| 2 . (20) 



The matrix S n is a n x n 
of row and column indices 



beplitz matrix (i.e. its matrix elements depend only on the difference 
lfl|). The generating function a(q) of a Toeplitz matrix is defined by 

(s„u = r % ^y q[3 ~ k) (2i) 







and in our case can be found from ()12I13|) as 

„(n\ = -(l + e"Q - - + — ° Sg ~ h + i7sing — 
2 ^ / 2 2a/(cos g - /i) 2 + 7 2 sin 2 g 

Thus, the problem of calculation of the EFP 

P{n) = |det(S n )| , (23) 

is reduced (exactly) to the calculation of the determinant of the n x n Toeplitz matrix S n defined 
by the generating function Q21I22[) . The representation (I23|) is exact and valid for any n. In our 
study we are interested in finding an asymptotic behavior of Q23|) at large n — > oo. 2 

Most of these results are derived using known theorems on the asymptotic behavior of Toeplitz 
determinants. We collect these theorems in Appendix^ In the following sections we apply them 
to extract the corresponding asymptotes of P{n) at n —* oo in the different regions of the phase 
diagram. Two major distinctions have to be made in this process. For the critical isotropic (7 = 0) 
XY model, one applies what is known as Widom's Theorem and one finds a Gaussian behavior 
with a power law prefactor [jjj . In the rest of the phase diagram, we apply different formulations of 
what is known in general as the Fisher-Hartwig conjecture, which always leads to an exponential 
asymptotic behavior. As expected, we find a pure exponential decay for the EFP in the non-critical 
regions. 

For h > 1, the exponential decay is modulated by an additional oscillatory behavior. 

At the critical magnetizations h = ±1, we discover an exponential decay with a power law 
pre-factor. Moreover, by extending the existing theorems on Toeplitz determinants beyond their 
range of applicability, for h = ±1 we obtain the first order corrections to the asymptotics as a 
faster decaying power law with the same exponential factor. For h = 1, the first order correction is 
also oscillating and this means that the EFP presents an oscillatory behavior of the EFP for h > 1. 

The reader who is not interested in the mathematical details of our derivations can find the 
results in Tableland skip the following sections to go directly to Sec. IVIIll where we analyze the 
crossover between the Gaussian behavior at 7 = and the asymptotic exponential decay at finite 
7 using a bosonization approach. 

2 The reader might notice that our generating function 1221 is almost the same as the one analyzed by Barouch 
et al. in [l^l (cri3] (q) = cos q ~ h+n slnq ). The only difference is the shift by the unity in our expression. 

(cos q — /l) 2 +"y 2 sin 2 q 

This difference changes dramatically the analytical structure of the generating function, in particular, its winding 
number around the origin, and requires a new analysis of the generated Toeplitz determinants. 
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(22) 



EFP for the Anisotropic XY model 


Region 


7, h 


P(n) 


Eq. 


Section 


Theorem 




h< -1 


fie-"* 3 




Iv All 


Szego 


So 


-l<h<l 






IV A 21 


FH 


£+ 


h > 1 


J5[l + (-l)M] e" n/3 




IV A 31 


gFH 


T E 


7 2 + h 2 = 1 






eh 


Exact 


n+ 


/i = 1 


£ n- 1 / 16 [1 + (-1)M/Vn] e" n/3 


Ell 


IVH1I 


gFH 




/i = -1 






IVB2I 


gFH 


n 


7 = 0, \h\ <1 


£;n- 1 / 4 e~" 2Q 




EH 


Widom 



TABLE I: Asymptotic behavior of the EFP in different regimes. The exponential decay rate j3 is given by 
Eq. lj25J) for all regimes. The regions in boldface are the critical ones. The coefficients E, A are functions 
of h and 7 whose explicit expressions are provided in the text. Relevant theorems on Toeplitz determinants 
are collected in thelAl 



(a): £_ (b): SI (c): £„ 




IV. SINGULARITIES OF a(q) AND EXPONENTIAL BEHAVIOR OF THE EFP 



To derive the asymptotic behavior of the EFP we rely on the theorems on determinants of 
Toeplitz matrices. These theorems depend greatly on the analytical structure of the generating 
function (|22|). especially on its zeros and singularities. 

Setting 7 = in (|22|). we see that for the Isotropic XY model the generating function has only 
a limited support within its period [0, 2ir}. This case is covered by what is known as Widom's 
Theorem and will be considered in Section fVHI 
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In the remaining parts of the phase-diagram the generating function has only pointwise singu- 
larities (zeros) as it is shown in Fig. [21 These cases are treated under a general (not yet completely 
proven) conjecture known as the Fisher-Hartwig conjecture (FH), which prescribes the leading 
asymptotic behavior of the Toeplitz determinant to be exponential in n: 



P(n) n ~°°e- pn . (24) 

While the pre-exponential factors depend upon the particulars of the singularities of the gener- 
ating function, the exponential decay rate is given in the whole phase diagram (7 7^ 0) according 
to FH as 



/?(/*, 7) 



f Z7r dq 

L 27 log|cx((?)| 

n dq 



2n l0g 




cos q — h 



(25) 



(cos q — h) 2 + 7 2 sin 2 q / 

The integral in (|25]> is convergent for all h and all 7 7^ and j3{h,^) is a continuous function of its 
parameters. 

In Fig. |31 (3(h, 7) is plotted as a function of h at several values of 7. One can see that 7) 
is continuous but has weak (logarithmic) singularities at h = ±1. This is one of the effects of the 
criticality of the model on the asymptotic behavior of EFP. 

These weak singularities are also a manifestation of the rich analytical structure underlying 
/3(/i,7) and the generating function (|22|) . To understand these structures, we interpret the periodic 
generating function (|22j) as the restriction to the unit circle (z = e 1 ®) of the complex function 

^) = |fl + ' ), (26) 

where 

Pi(z) = ^(z- Zl )(z-z 2 ), (27) 
P2(z) = ^L(z lZ -l)(z 2 z-l) (28) 

with 



22 = W^2Zi. (30) 

The integral in (|25|) can be regarded as a contour integral over the unit circle of the function 
(|26|) . We can deform the contour of integration taking into account the complex structure of the 




FIG. 3: Plot of the decay rate f3 as a function of the parameters 7 and h. The function diverges for 7 = 
and is continuous for h = ±1 (although has weak singularities at h = ±1). 

integrand in the various regions (see Fig. 0J) and express (|25j) as a simpler integral on the real axis 
(after partial integration and some algebra). 



A. The non-critical regions (S± and S ) 



1. £_ (h < -I) 



For h < —1, the analytical structure of the integrand of (|25|) is shown in Fig. 0^,. We re-write 
the decay rate ()25|) in this region as 



7) = 2 ln 



yjh 2 + 7 2 - 1 - h 
7 + 1 



A^jj-A^T), 



(31) 



where 



A(fc, 7 ) = In 
A(/i,7) = 



1 



1 



/i /l -7 



2 V' |/i|Vl + 7 
1 dx 1 



\K\ 2vr ^(l-x 2 )^ 2 -^ 2 ) V z 



x H m 



x — a 



x + a 



(32) 
(33) 
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FIG. 4: The integral in (|25|) is performed over the unit circle C. The analytical structure of the integrand 
allows for a deformation of the contour of integration into T, which encloses a logarithmic branching line, 
different in the various regions of the phase-diagram (in (d), T encloses also a simple pole at the origin). 
The roots z\ and zi were defined in ^ and Ipp . 



with 



K = 



yjh? + 7 2 - 1 - 7 

y^ 2 + 7 2 - 1 - 7 
/i-l 



(34) 
(35) 



This decomposition of /3(h, 7) is especially useful in analyzing the transitions between non-critical 
and critical regimes. In fact, we will see that the functions A(h, 7) and A(/t, 7) denned above are 
universal across the phase diagram (hence the need for the seemingly redundant absolute values in 
our definitions). 



2. So (\h\ < I) 

As before, we can express the contour integral defining /3(h, 7) as a standard integral on the 
real axis. For \h\ < 1 and h? +7 2 > 1, the structure of the integrand is depicted in Fig. and the 
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decay rate is simply 

/3(/», 7 ) = -A(/i j7 )-A(/»,7), (36) 

where A(h, 7) and A(/i, 7) have already been defined in (|32I33|) . 

For h 2 + 7 2 < 1, the structure is quite different (see Fig. 011). In this region the expression for 
/3(h, 7) in terms of a real axis integral is complicated and will therefore be omitted it in this paper. 



3. £+ (h > 1) 

A calculation similar to the previous ones (see Fig. 0f ) gives the expression for the decay factor 
for h > 1: 

/?(/>, 7) = ^ In 

where A(h, 7) and A(fo, 7) were introduced in (f3*2")l and (|3"3*|) . 

One important difference exists in this region: as will be discussed in length later in Section 
IV A 31 in S + there are two equivalent representations of the generating function. This ambiguity 
reflects on the value of (3, in that the choice of the representation for the generating function 
determines the branch cuts in Fig. 0J We will see that we have to use both values of (3, which differ 
only by an imaginary constant: 

(3' = f3 + m (38) 
and this will add an oscillatory behavior to the EFP. 



7 + 1 



1 + h 



A(fc,7)-A(/i, 7 ), 



(37) 



B. The critical lines (ft±) 



We can calculate the decay factor f3 at h = 1 (fi+) from a limiting procedure on ()36|) or (|37|). 
At h = 1, only A(/i, 7) is nonvanishing, thus guaranteeing the continuity of (3 across the critical 
line. From an appropriate limit of (|33|) . we calculate the decay rate for h = 1 as 

~ l dx 1 



0(1,7) 



In 



1 — jx 



hi 



1-7 
1 + 7 



(39) 



1 + 73; 

As discussed before in reference to (|3*7)l . the definition of 7) is not unique and, as in l|38jl. will 
generate again an oscillatory behavior for the EFP (see later in Sec. IV B"T|) . 

The value of /3 at h = — 1 can also be obtained from a limiting procedure on (|3"3"|) 

" x dx 1 



/?(-!, 7) 



2vr^T 



:hl 



1 — 7X 



1 + 7X 



In 



1 + 



1-7 
1 + 7 



(40) 
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As can be seen in Fig. |HJ the decay factor f3 is continuous across the critical lines, but has a 
discontinuity in its derivative. As (5 approaches the critical lines, it actually shows a non-analytical 
behavior leading to a logarithmic singularity: 



The derivative d(3/dh diverges logarithmically as h — > dbl. 

Moreover, one can easily notice from the difference between expression (|36|) and 1)3 7 Jl that even 
the finite part of the derivative of /3(h, 7) by h is different if one approaches the critical line h = 1 
from above or below, due to the appearance of the additional term in 1)37|) . The same holds across 
the critical line h = — 1, due to the presence of the first term in ()31j) . which doesn't appear in 1)36)1 . 

V. THE PRE-EXPONENTIAL FACTORS 

For 7 ^ 0, the leading behavior of the EFP is always exponential. However, the singularities of 
the generating function are different in different regions of the phase diagram and we must therefore 
use different forms of the Fisher-Hartwig conjecture in order to derive the pre-exponential factors 
and determine the asymptotic behavior of P{n). We will now show how to obtain the results for 
each of the regions. 

A. The non-critical regions (S± and S ) 

1. £_ (h <-l) 

In this region (7 ^ 0, h < —1) the generating function 1)22)) is nonzero for all q (see Fig.|2ti): this 
is the simplest case and can be treated using the (rigorously proven) Strong Szego Limit Theorem, 
see (|A3|) . It gives 



P{h = ±1 + 6,7) = /?(±1,7) + - eln|e 



(41) 



P{n) = |det(S n )| 



OG 



in 



(42) 



with /3(/i, 7) given by (|23)) or (|3"T|) and 




(43) 



where is defined in (|A5|) as the fe-th Fourier component of the logarithm of a: 
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° V V ( cos 1 ~ ^) + 7 2 9 / 

The sum in (|43|) is convergent only for 7 7^ and for h < — 1. For /i > — 1, the generating function 

(I22|) develops singularities which produce 1/k contributions to ()44l) that make the sum in (J43|) 

divergent. Therefore, in the rest of the phase diagram these singularities have to be treated to 

absorb the harmonic series contributions. Consequently, each region of the phase diagram will 

involve a different definition for the pre-exponential factor and the " regularization" procedure will 

sometimes generate an additional power-law contribution. The result is given by the Fisher-Hartwig 

conjecture that we must use in the remainder of the phase diagram. 

2. S (\h\ < 1) 

As can be noticed from Fig.[2t, in So (7 7^ 0, —1 < ft, < 1) the generating function a(q) vanishes 
and its phase has a discontinuity of tt at q = tt. The asymptotic behavior of Toeplitz determinants 
with this type of singularities in the generating function is given by FH, which is actually proven 
for cases in which only one singularity is present. 

We decompose the generating function as in (|A6|) 

a(q) = T{q)M q -^ mod 2n ~^ (2-2 cos(g - ir)) 1/2 (45) 
and using (|A7|) we obtain 

P{n) = |det(S n )| n ~°°Eo(h, 7 ) e ~^ n . (46) 

The behavior is exponential as before with the decay rate (3{h,~f) from (|25I36|) . but the pre- 
exponential factor is different. According to (|A8|) it is given by 

Eo (h, 7 ) = ^- (47) 

where, as in (|A4j) and (|A5|) 

E[t] = exp (jrkT k T-k) (48) 



\k=i 



and 



(-l) k 

f k = a k - 6(k). (49) 

Here 9{k) is the usual Heaviside step function. As we mentioned in the previous section, a k (|44|) 
has 1/k contributions from singularities of cr(q) and the effect of the parametrization ()45l) is to 
cure (remove) these harmonic series divergences of the prefactor of the EFP in this regime. 
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3. £+ (h > I) 

In E_|_ (7 ^ 0, h > 1), <r((7) vanishes at q = and q = ir and its phase presents two 7r jumps at 
those points (Fig. 

In this case the application of FH leads to some ambiguity, because there exist two representa- 
tions of the kind (|A6|) and one obtains two values for /3(h, 7) using the two representations of the 
generating function: j3\ = j3 and P2 = P + ivr, with [5 from ((251) or (|37j) . This ambiguity is resolved 
by the (yet unproven) generalized Fisher-Hartwig conjecture (gFH), which gives EFP as a sum of 
two terms so that both values of /3's are used (see the appendix I A 31 or Q]). 

The two leading inequivalent parametrizations ((A13|) are: 

a (q) = r 1 (g)e5[(^) mod ^- 7r ](2-2cos(g-^)) 1/2 

xe -i[,mod2^] ( 2 -2cosq) 1/2 (50) 

= T 2 {q)e-^ q ~^ mod (2-2 cos(g - vr)) 1/2 

Xe i[?mod2 7 r- 7 r]( 2 _ 2cosg > ) l/2 ^ 

Application of (|A14I) gives the asymptotic behavior of the determinants as 

|det(S n )f~°° [El(h n ) + (-l)» E *(h,>y)] e~^ n (52) 



with 



E*{h,i)= f [T] (53) 
t+(0)t_(tt) 



^IC.7)-^fL (54) 
and /3(/i,7), £[r] defined in (I23HH1) with 

r fc = a k - t^L e{k) - i (55) 

Once again, as in the previous section, the effect of the parametrization is to remove the 1/k 
contributions to (|44|1 due to the singularities of the generating function. 

We conclude that the non-critical theory presents an exponential asymptotic behavior of the 
EFP. In the region E + , however, the EFP in addition has even-odd oscillations 

P{n) n ~°°Ei{h, 7) [1 + A+(h, 7) cos(vrn)] e"^'^, (56) 

where the exponential decay factor is given by (|37|) . 
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The amplitude of the oscillations is 



t_(0)t + (tt) 
t(0) /t_(tt) x 2 



A + (h,j) 



r(vr) V r -(0) 

( av I dz lo S T ( z ) \ , K -s 
= —I 6XP { 4 l™of to z*-(l + e)* ) ' (57) 

where we used (|A9|) . the definition of r and ()A11|) . We can deform the contour of integration as 

in Fig. Ilf and calculate the integral in Q57JI to obtain 

, Jh 2 - 1 

A+(/i,7) = 0r(Vy)= 7 , 2 2 = — , (58) 

where if (ft,, 7) was defined in (|34[). 

Expression (|56|) for the EFP fits the numerical data remarkably well (see Fig. |SJ) and this fact 
strongly supports the generalized Fisher-Hartwig conjecture. 

One can understand these oscillations as a result of "superconducting" correlations of real 
fermions described by the Hamiltonian (J7J). Fermions are created and destroyed in pairs of nearest 
neighbors. At large magnetic fields, the oscillations are due to the fact that the probability of 
having a depletion string of length 2k — 1 or 2k is very similar. Since the magnetic field in is 
essentially a chemical potential for the fermions, the energy cost to destroy a pair of particles is 
Ah: at very big magnetic fields, the amplitude for a pair destruction event is suppressed by a factor 

2 

of i.e. a probability of jgp- This means that the probability of depletion behaves like: 

/4h\~ 2k 
P(2fc-1) ~ 2 I —j and 

P{2k) ~ ^ ^ , (59) 

where the factor of two in the first expression is a simple combinatorial factor. The two probabilities 
in (|59|) can be combined in a single expression: 

—J , (60) 



1-T + ofA). (Oil 



which is precisely (|56j). with 

We can check the correctness of this interpretation by taking the limit of (|56l) for h >> 1,7. 
From and l(5S|) it is easy to find 

^-,00,7) = log^+ (62) 
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FIG. 5: Results of the numeric calculation of the Toeplitz determinant are shown as P(n)e^ n vs. n at 
7 = 1, h = 1.1. The value of f3 is obtained numerically from l|39fl . The solid line is the analytic result 
E{\ + (-l)'M) with A = 0.2182... from J5J) and E = 0.6659... obtained by fitting at large n. To make the 
plot more readable we show only every 11th point (for n = 1, 12, 23, . . .) of the numerical calculation of the 
determinant. Note that the size of the points is not related to the estimated error in the numerics, which is 
actually smaller. 



A + (h 00,7) 



(63) 



in agreement with ()6UI61|) . 



B. The critical lines (£1±) 

1. n+ (h = 1) 

For h = 1 the generating function cr(q) vanishes at q = it and its phase has tt jumps at q = 0, n 
(see Fig. HJl). As in the previous section, the existence of two singular points gives rise to many 
terms of the form (|A13|) , However, in contrast to the S + region, the application of gFH as in (|A14|) 
shows that all terms are suppressed by power law factors of n with respect to the leading one. 
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The leading term is generated by the parametrization: 

a{q) = T 1 (q)e^ q -^ mod27T -^ (2 - 2cos(q - tt)) 1/2 e-^ gmod2 ^ (64) 

and consists of an exponential decay with /?(1, 7 ) from (|39j) and a power law contribution with 
critical exponent A = 

|det(S n )| ~ E\ (j)n-TS e -M,7> (65) 

with 

E\ ( 7 ) = E[t]G ( 7 Wf) , (66) 

L W W 2i/4 r i/ 4 ( )r_(7T) 

where G is the Barnes G- function defined in (|A10|) and E'fr] is defined as in (1381) with 

r fe = a fc + Q - (-l) fc ) 1 - ^ 9{—k), (67) 

with <Tfc from (|44|). 

The next term (subleading at n -> oo) is obtained from the parametrization 

a{q) = T\q)e-^ q - w)mod2w - w] (2 - 2cos(q - tt)) 1/2 e^ [qmod27T - n] (68) 

and is given by 

Ef {l){-l) n n~TE e-PQ-ifr ( 69 ) 

with 

B!w S Ew(i)aa) f m ■ (70) 

Although the inclusion of the latter (subleading) term is somewhat beyond even gFH, we write 
the sum of these two terms as a conjecture for EFP at h = 1 

P(n) ~ E\ (7) U + (-l) n A 1 ( 7 )/n^j e"^ 1 ' 7 ^. (71) 

As these results rely on our unproven conjecture, we present our numerical data for this case 
in Fig. 113 Indeed, we see that the form (|71|) is in good agreement with numerics and this supports 
our hypothesis. 

The amplitude of the oscillations is 

Al[V ~ V2 G(1)G(|) r_(0)r + (vr) 
1 r(f) r(0) /r_(7r)' " 



72 T(i) r(vr) Vr_(0) 
£itl I exp (4 lim / d " lQgT ^) ^ (72) 
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FIG. 6: Results of the numeric calculation of the Toeplitz determinant are shown as P(n)e' 3n n 1 / 16 vs. n 
at 7 = 1, h = 1. The value (3 = log 2 + 2G/tt with Catalan's constant G is obtained from (|25[l . The solid 
line is the analytic result £7(1 + (-l) n A/n?) with A = 0.2399... from O and E = 0.8065... as obtained by 
fitting at large n. To make the plot more readable we show only every 11th point (for n = 1, 12, 23, . . .) of 
the numerical results on the determinant. Note that the size of the points is not related to the estimated 
error in the numerics, which is actually smaller. 



where we used (|A9|) and the identity 



G(z + 1) =T(z)G(z) 



(73) 



To calculate the integral we deform the contour of integration as in Fig. and find 

1 



^4l(7) 



£(i: 
r(j; 



(74) 

oo but with an additional 



We conclude that at h = 1 the EFP decays exponentially at n 
power law pre-factor and a damped oscillatory component. 

Remark. It is curious to notice that the exponents 1/16 and 9/16 in (|65|) and 1)69(1 remind us 
of the scaling dimensions of spins a x and a y . 3 It looks as if the EFP operator @, among other 
things, has inserted square roots of spins transverse to the magnetic field at the ends of the string. 



3 See Ref. |l3j|. where it was shown that the power laws for the a x and a y correlators are 1/4 and 9/4 respectively. 
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2. 0_ (h = -l) 



For h = — 1 the generating function u{q) does not vanish but has a phase discontinuity of n at 
q = tt. We parametrize cr(g) as 



and apply FH to obtain 



with 



a(q) = T 1 (q)e-^ q ' 7T)mod2n - n] 



P{n) ~ £^(7) e-^f- 1 '^" 



00) 



where /3(— 1,7) and E[r] are defined in (|4()|1 and (|48j) with 



(75) 



(76) 



(77) 



(78) 



4fe 4fc 
and from (|44|l. 

We can stretch the gFH the same way as in the previous section for h = +1 by considering the 
second parametrization 



a 



(q\ = r 2 (g)e i ^ 9 ~ 7r ) mod27r ~ 7r ] 



which gives 



with 



P'(n) ~ £^(7) n"T6 e" 



(7) = E[t]G 7 C 7 



-/3(-l, 7 )n 



4J ~ \AJ r 3/4 (?r)) 

Adding this subleading term to (|76jl we obtain 

P(n)~^i J (7)n - ^ fl + A_i( 7 )/ra^ P-^-L-r) 



(79) 



(80) 



(81) 



(82) 



with 



A-i( 7 ) 



g(|)g(D r + (7r) 
r(|) r + (vr) 



r(|) r_(7r)- 

We propose (|82|) as an asymptotic form for EFP at h = — 1. 



(83) 
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VI. THE LINE T E : AN EXACT CALCULATION 



Before we conclude our analysis of the EFP with the study of the isotropic XY model, let us 
check our results Q42I25|) on the special line 4 in the phase diagram defined by 

/i 2 + 7 2 = 1. (84) 

It was shown in Ref. ^1 that on this line the ground state is a product of single spin states and 
is given by 



ig> = n = n [ c ° s (}) i t>i) + (0 1 1, j) 



(85) 



where | |, j) is an up-spin state at the lattice site j, etc. One can directly check that the state (|53|) 
is an eigenstate of @ if the value of parameter 6 is 

cos 2 9 = (86) 

n 1+7 

and ()84j) is satisfied. It is also easy to show [Tj] that this state is, in fact, the ground state of Q. 
The probability of formation of a ferromagnetic string in the state (|85|) is obviously 



which is an exact result on the line ()84l) . The value of /3{h, 7) which immediately follows from this 
exact result is 

[3(h = ±^T^, 7 ) = - log Q T \\Jt^) = ~Mh,l), (88) 

where A (h,j) was defined in (|32|) . 

This is, indeed, consistent with ()36|) since under the condition ()84j) the function A(/i, 7) vanishes. 
The integral (|33|) defining A(/i, 7 ) vanishes for ()84() because the branching points ()29|) and (|3(J|) 
collapse to the same point and therefore the region of integration shrinks to just one point (|33|) . In 
fact, the Toeplitz matrix (|21j) generated by (1221) becomes triangular on the line (|84|) with diagonal 
matrix element (S n )jj = sin 2 (0/2) and the determinant of S n is exactly (j87|) . 

From the definitions of /3(/i,7), we see that the decay factor consists of two terms, which have 
now a clear physical meaning. The A(/i, 7) term is the factor we found above in (|88|) and represents 
the contribution given by un-entangled spins. The remaining part accounts for the correlations 
between spins. Both A(/i, 7) and the correlation functions given by (|12|) and ()13|) vanish on the 
line (jHU. 



4 We are grateful to Fabian Essler who suggested us to check our results on this special line and pointed out the 



reference 01 to us. 
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VII. THE CRITICAL LINE n (7 = 0) AND THE GAUSSIAN BEHAVIOR 

The case 7 = 0, corresponding to the Isotropic XY Model, has been studied in Ref. Q|. For 
7 = the generating function ()22j) is reduced to the one found in Q]. 

For \h\ < 1, the generating function a(q) has a limited support between [— cos -1 h, cos -1 h]. To 
find the asymptotic behavior of the determinant of the Toeplitz matrix one can apply Widom's 
Theorem 2!) and obtain [?| 



2 



P{n) ~ 2Ae 3 «'(- 1 )(l - h)-ln-\ ' " (89) 

We see that in this case, the EFP decays as a Gaussian with an additional power-law pre-factor. 

In a different context, the formula (|89j) appeared also in probability of forming a gap 

in the spectrum of unitary random matrices. This is not unexpected, since the joint eigenvalue 
distribution of unitary random matrices is known to coincide with the distribution of free fermions 
in the ground state. 

For \h\ > 1, the theory is no longer critical and the ground state is completely polarized in the 
z direction, giving a trivial EFP P{n) = for h > 1 and P(n) = 1 for h < — 1. 



VIII. CROSSOVER BETWEEN GAUSSIAN AND EXPONENTIAL BEHAVIOR: A 
BOSONIZATION APPROACH 

In order to understand qualitatively the crossover between the Gaussian asymptotic behavior 
at 7 = and the exponential decay for 7 ^ 0, we employ a bosonization approach similar to the 
one used in In the limit 7 <C 1 we consider the continuum limit of (JHJ), bosonize the fermionic 
fields, and write the Euclidean action of the theory as S = J dx drC, where r = it is the imaginary 
time and the Lagrangian is 



C = 2 sin kp 



(d^ f - - cos (Vi7n?) • (90) 



7T 

Here hp = cos -1 h is the Fermi momentum at 7 = 0. 

This is a Sine-Gordon theory for the "conjugate field" -!9(x,r), which describes the imaginary 
time dynamics of our 1-D system. In terms of $ the density of fermions is given by p = d T i!) + po, 
where po = kp/n is the density of fermions in the ground state. 

In the field theory approach, the EFP (see Ref. (|) in the limit n —> 00 would be given with 
exponential accuracy by the probability of an instanton P(n) ~ e~ So , where So is the action of 
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the instanton. Here the instanton is the solution of the classical equations of motion of (|9U|) which 
corresponds to the formation of an emptiness of length n at the time r = 0. Unfortunately, the 
EFP instanton involves large deviations of the density of fermions from the equilibrium density po 
and is beyond the bosonization approach as the derivation of (j9U|) relies on the linearization of the 
fermionic spectrum near the Fermi points. 

Following [6||, we are going to slightly generalize our problem, by considering the depletion 
formation probability instead of the EFP requiring 

P\t=0, 0<x<n = P0+ df&fa t)\ t =0, 0<x<n = P0 ~ P, (91) 

where p is some constant. The original EFP problem corresponds to p = pq. Here, instead, we 
consider the probability of weak depletion, i.e. 

- d t #(x,t)\ t=0i0<x<n = p«p . (92) 

We study the latter using an instanton approach to ()9U|) and infer the (qualitative) behavior of the 
original EFP from this weak limit. 

To simplify the problem further, we assume that the instanton configuration is completely 
confined to one of the wells of the Cosine potential in (|9U|) and that the field i? is small enough to 
allow for an expansion of the Cosine: 

S ~ 2 sin kF I dx dr 



J dx dr (d^f + 2 7 tf 2 . (93) 



In this formulation, the anisotropy parameter 7 1 / 2 plays the role of the mass of the bosonized theory. 
The probability we are looking for is given by the action Sq of the classical field configuration which 
satisfies the Euler-Lagrange equation (in this case a Klein-Gordon equation in two dimensions) with 
the boundary condition (|9~Tj) 

Pp{n)=e- S ". (94) 



In the limit 7 = 0, the theory is massless and scale invariant. In |6j it was shown that, due to 
the scale invariance, the action of the instanton is quadratic in n. The instanton configuration in 
this case is essentially a droplet of depletion in space-time with dimensions proportional to n both 
in the space and time direction, in order to satisfy the boundary condition Q91|) . This result is 
consistent with the Gaussian asymptotic behavior prescribed by Widom's theorem (see Sec. IVlljl . 

In the massive case, a finite correlation length £ ~ ^~ 1 / 2 is generated and one has a crossover 
behavior. For string lengths n smaller than the correlation length 7 _1//2 , the instanton action is 
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FIG. 7: Plot of the value of the stationary action Sq FIG. 8: Plot of the derivative dSo/dn with Sq from 



vs. the string length n. The action Sq is obtained (IG9|) . The plot corresponds to to = 0.01, p = 0.2 and 
from (|G9p with f(y) given by the numerical solution clearly shows a crossover from the quadratic to the 
of the singular integral equation l)C7|) . The graph linear behavior at n ~ 2/to = 200. 



depicts S (n) for to = ^ = 0.01, p = 0.2. The 
crossover takes place around n ~ 2/m = 200. 



not sensitive to the presence of the finite correlation length and is still quadratic in n (giving a 
Gaussian decay for EFP). In the asymptotic limit of string lengths greater than 7 -1 / 2 , the time 
dimension of a depletion droplet is of the order of £ (instead of n as in the massless limit): the 
action is linear in n and the probability has an exponential behavior. 5 

In Appendix [C] we show how to solve the integral equation corresponding to the boundary 
problem (|91I93[) and present its numerical solution and some analytical results. Figures and |H1 
clearly show the crossover between a quadratic behavior of the stationary action for small n to a 
linear asymptotic one for n — > 00. 



IX. DISCUSSION AND CONCLUSIONS 



The asymptotic behavior of the Emptiness Formation Probability P(n) as n — > 00 for the 
Anisotropic XY model in a transverse magnetic field as a function of the anisotropy 7 and the 
magnetic field h has been studied. We have summarized our results in Table [IJ These asymptotic 
behaviors have already been presented in ^J. In this work, we completed the derivations by 
providing explicit expressions for the coefficients of these asymptotic forms. 

Our main motivation has been to study the relation between the criticality of the theory and 

5 This picture is very similar to the one for massless theory at finite temperature. In the latter the inverse temperature 
plays the role of the correlation length Q fsee lc)l . 



24 



the asymptotics of the EFP. Let us now consider the results on the critical lines (Oo and fi-jj. The 
Gaussian behavior on SIq (7 = 0, \h\ < 1) is in accord with the qualitative argument of Ref. [6j using 
a field theory approach. In So (7 7^ 0, \h\ < 1) the asymptotic decay is exponential. We proposed a 
physical interpretation of the crossover between the two asymptotes using a bosonization analysis 
of the region of small 7: we suggest that there is an intermediate regime of Gaussian decay for the 
string lengths smaller than 1/^/7 which crosses over to the exponential behavior for longer strings. 

On the critical lines Q±, the decay of the EFP is exponential instead of Gaussian, and apparently 
contradicts the qualitative picture of Ref. jg|. The reason for this disagreement is that although 
at h = ±1 the model can be rewritten in terms of massless quasiparticles x defined in Q, we 
are still interested in the EFP for the "original" Jordan- Wigner fermions ip. In terms of x this 
correlator has a complicated (nonlocal) expression very much different from the simple one (|14j) . 
From the technical point of view, the difference is that in the qualitative argument in favor of 
a Gaussian decay of EFP for critical systems there is an implicit assumption that the density of 
fermions (or magnetization) is related in a local way to the field responsible for the critical degrees 
of freedom (free boson field <p). This assumption is not valid on the lines h = ±1. The theory is 
critical on those lines and can be described by some free field <j>. However, the relation between 
the magnetization and this field is highly nonlocal and one can not apply the simple argument of 
f| to the XY model at h = ±1. 

Although EFP at the critical magnetic field does not show a Gaussian behavior, there is an 
important difference between the asymptotic behavior of EFP on and off critical lines. Namely, a 
power-law pre-factor n~ x appears on all critical lines. For the XY model it is universal (i.e. A is 
constant on a given critical line) and takes values A = 1/4 for 7 = :7] and A = 1/16 on the lines 
h = ±1. It would be interesting to understand which operators determine these particular "scaling 
dimensions" of the EFP (see the remark at the end of Section IV B 1|) . 

At h > 1 the use of gFH predicts even-odd oscillations of P(n). We compared the predicted 
oscillations to numerical calculations of Toeplitz determinants and found a very good agreement 
(see Figs. I">l(i|) . We proposed a physical interpretation of the oscillations as coming from pair 
correlations of spins which can be clearly seen as superconducting correlations in the fermionic 
representation ((7J). 

In some parts of the phase diagram (S + , £1±) we used the so-called generalized Fisher- Hartwig 
conjecture Ijj] which is not yet proven. However, our numeric calculations support the analytical 
results (see Figures 03 and H)J) . We note that to the best of our knowledge this is the first physically 
motivated example where the original Fisher-Hartwig conjecture fails and its extended version is 
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necessary. 6 We also suggest that the gFH could be used to find the subleading corrections to the 
asymptotic behavior, as we did for h = ±1 in (|71I82|) . This novel hypothesis is supported by our 
numerics and it would be interesting to confirm it analytically 

In conclusion, we notice that it is straightforward to generalize our results for nonzero tem- 
perature. The only modification is that at T / the thermal correlation functions must be used 
instead of (|12I13() . Then, the generating function (|22|) is non-singular everywhere and we have an 
exponential decay of P(n) in the whole phase diagram according to the standard Szego Theorem 
and standard statistical mechanics arguments. We present results for T ^ in Appendix iBl 
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APPENDIX A: ASYMPTOTIC BEHAVIOR OF TOEPLITZ DETERMINANTS 

The asymptotic behavior of the EFP for at n — > oo is exactly related to the asymptotic 
behavior of the determinant of the corresponding Toeplitz matrix (|21|22|23|) and can be extracted 
from known theorems and conjectures in the theory of Toeplitz matrices. These types of calculations 
have been done first in ["Lsl . Iisl | for spin-spin correlation functions. It is well known that the 
asymptotic behavior of the determinant of a Toeplitz matrix as the size of the matrix goes to 
infinity strongly depends upon the zeros and singularities of the generating function of the matrix. 



23j. Here we 



A very good report on the subject has been recently compiled by T. Ehrhardt 
want to recap what is known about the determinant of a Toeplitz matrix generated by a function 
a(q): 

D n [a] = det(S n ) = det ( f a^e^^p-Y , (Al) 
where the generating function a(q) is a periodic (complex) function, i.e. cr(q) = <t(27t + q). In this 



6 We note that recently the theory on Toeplitz determinants has been used and extended with new results in order 
to calculate yet one more important physical quantity. We refer the interested reader to [2tT |. [ill ] and |2^ |. where 
the entaglement for the XY Spin chain and for Random matrix models have been calculated. 
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work we dealt only with generating functions with zero winding number 

lnda(q)^ £ ^ A l og(T (g) = (A2) 

and this brief review will be limited to this condition. This was not the case in the study of Barouch 
et al. where the generating function (see footnote after (j2HJ)) had non-zero winding number 
in some regions of the phase diagram. 



1. The Strong Szego Theorem 

If a(q) is sufficiently smooth, non-zero and satisfies Ind cr(q) = (he-, the winding number is 



251] h which states that 



0), we can apply what is known as the Strong Szego Limit Theorem ( 
the determinant has a simple exponential asymptotic form 

D n [a] ~ E[a]G[a] n n -» oo, (A3) 

where G\a\ and E[a] are defined by 

oo 

G[a] = exp o"o, E[a] = exp k&ko~-k (A4) 

k=l 

and are the Fourier coefficients of the expansion of the logarithm of o~(q): 

oo 

log a(q) = *ke ikq . (A5) 

k=— oo 

2. The Fisher-Hartwig Conjecture 

Over the years, the Szego Theorem has been extended to consider broader classes of generating 
functions by relaxing the continuity conditions which define a " smooth function" , but it remained 
limited to never-vanishing functions. Therefore, some extensions have been proposed to the Szego 
Theorem in order to relax this latter hypothesis. When the generating function has only pointwise 
sing ularities (or zeros), there exists a conjecture known as the Fisher-Hartwig Conjecture (FH) 

7 

When a(q) has R singularities at q = 6 r (r = 1..R), we decompose it as follows: 

R 

a(q) = T {q) \\ e ^[(<i-e r ) mod a™] (2 _ 2 cog(g _ q^K (A6) 



This conjecture is still not completely proven. For details and status of the conjecture see Ref. [ig|. 
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so that r(q) is a smooth function satisfying the conditions stated in the previous section. Then 
according to FH the asymptotic formula for the determinant takes the form 

D n [a] ~E[T,{K a },{*ah{6a}] n^-l) G[r] n n ^ 0O) (A7) 

where the constant prefactor is conjectured to be 



r=l 

(« r +Ar)(«s— A s ) 



l<r^s<R 

^ G(l + K r + \ r )G(l-K r + \ r ) 

1J WTv^) • (A8) 

E[t] and G[t] are defined as in (|A4I) and r± are defined by decomposition 

r(g)=T„(e^)G[r]r + (e-^), (A9) 

so that r + (r_) are analytic and non-zero inside (outside) the unit circle on which r is defined 
and satisfy the boundary conditions t+(0) = r_(oo) = 1. G is the Barnes G-function, an analytic 
entire function defined as 

oo , 
,21 /o TT /„ ■ - 2 

n 



of a single singularity (R = 1), see 
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G(* + 1) = ( 2vr )-/2 e -[-+(7 B +i)^/2 H{l + ±\ e~ z +^ , (AlO) 

n=l n 

where "Je ~ 0.57721 ... is the Euler-Mascheroni Constant. 

This conjecture is actually proven for some ranges of parameters K r and X r or fully for the case 

In many simple cases it is possible to find the factorization of r into the product of t + and r_ 

by inspection. More complicated examples like the ones presented in this work require a special 

technique to obtain this factorization, which is known as the Wiener-Hopf decomposition: 

f dz logr(z) 
J 2tti z — w 

logr_H=-/— H>1, (All) 

J 27T1 Z — W 

where the integral is taken over the unit circle. 

In light of these formulas, it is useful to present the parametrization (|A6|) in a form which makes 
the analytical structure more apparent. Changing the variable dependence from q to z = e 19 , we 
write 

°(z)=r(z)H(l-^) (l-^) , (A12) 

where z r = e l9r . 



r=l v 
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3. The Generalized Fisher-Hartwig Conjecture 



Despite the considerable success of the Fisher-Hartwig Conjecture, few examples have been 
reported in the mathematical literature that do not fit this result. These examples share the 
characteristics that inequivalent representations of the form (|A6|) exist for the generating function 
cr(q). Although no theorem has been proven concerning these generalization of the Fisher- 

Hartwig Conjecture (gFH) has been suggested by Basor and Tracy that has no counter- 
examples yet. 

If more than one parametrization of the kind (|A6|) exists, we write them all as 

R 

a(q) = T \q) J] e iK ^~ 6 ^ mod 2 ^ (2-2 cos(g - 6 r )) K , (A13) 

r=l 

where the index i labels different parametrizations (for R > 1 there can be only a countable number 
of different parametrizations of this kind). Then the asymptotic formula for the determinant is 

A^M-X^trSRMAiMM] n^G[rT n - oo, (A14) 



where 



R 

m = E(( A ') 2 -(4) 2 ), ( A15 ) 



r=l 



T = U || Re [fi(»)] = maxRe [«(?)] j . (A16) 

The generalization essentially gives the asymptotics of the Toeplitz determinant as a sum of 
(FH) asymptotics calculated separately for different leading (see Eq. (|A16|0 representations l)A13|) . 
In Sec. IV BH we used the sum of all (not necessarily leading) representations and showed that it also 
correctly produces the first subleading corrections to the asymptotics of our Toeplitz determinant. 

4. Widom's Theorem 

If a(q) is supported only in the interval a < q < 2tt — a as in our model for 7 = 0, singularities 
are no longer pointwise and one should apply Widom's Theorem [gjj. It states that the asymptotic 
behavior of the determinant in this case is 

D n [a] ~ 2 1 /^ e 3C'(-D ^y 1 ' 4 E[ P fn-^G[ P T [co^f , (A17) 

where E and G are defined in (|A4|) and 

a 



p{q) = a \2 cos 1 cos — cos g J (A18) 
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with the convention < cos -1 x < tt. 

For the case considered in Section IVlll the generating function is constant, E[p] = G[p] = 1, 
and QA17JI simplifies considerably giving 

D n [a] ~ 2 1 /i2 e 3C'(-D ( sin |)" 1/4 n -V4 (cos I)"' . (A19) 

APPENDIX B: EMPTINESS FORMATION PROBABILITY AT FINITE TEMPERA- 
TURE 

At finite temperature (T > 0), the correlators (|12j) and (|13j) become 

Ff k = i(^>r = -i(VW>r = ^^m^tanh^e^^), (Bl) 

" 

•/ 

The EFP is expressed by © and in the spinless fermion formalism it becomes 

n 

P(n) = (H^4) T . (B3) 



i=l 



We again use Wick s Theorem (or its thermal version, called Bloch-de Dominicis theorem |3Q|) 
to express it as a Pfaffian. The calculation proceeds the same way as for zero temperature and the 
EFP can be represented as 

P{n) = |det(T n )|, (B4) 
where T n is the n x n Toeplitz matrix generated by the function 

t(q) = I Ci + e^Hanh^) (B5) 



2 V 2T 

where the "rotation angle" i? 9 and the spectrum e q were defined in (|l(Jj) and f| 1 If) respectively. 

The generating function t(q) is never-vanishing and has zero winding number. Therefore, for 
T > we can apply the standard Szego Theorem to obtain 

P(n) n ~ 00 J B(/i,7,T)e- n/3(/l ' 7 ' T) , (B6) 

where 

r-'lK 



P(h,7,T) =~JJ^ log |*(^)| 
1 f 2 * dq 



— log 
2 7 2?r 6 



l(l+ COSg -\anh^ 
2 V e q 2T 



(B7) 

E(h, 7, T) = exp ( ki ki~k ) (B8) 



vfc=i 
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with 



, 2vr 

and e q is given as in (fTT|) by 



If cos q - h + ry sin q e q 

- 1 H tanh — =- 

£„ 2T 



e q = a/ (cos g — h) 2 + 7 2 sin 2 q. 



(B9) 



(BIO) 



As can be expected from simple thermodynamic considerations, at finite temperature the be- 
havior is always purely exponential asymptotically. As it was shown in [6j], at finite but very 
low temperatures one can observe a crossover from the zero temperature behavior at short string 
lengths n to the exponential behavior (|B6|) in the limit of very large n. This crossover occurs at a 
length scale of the order of the inverse temperature. 



APPENDIX C: CALCULATION OF THE STATIONARY ACTION IN THE BOSONIZA- 
TION APPROACH 

In Section IVIIII we have formulated the XY model near 7 = in terms of the bosonic field 
with Lagrangian ()93|) , It was also pointed out that, instead of the EFP, we are interested in 
the Probability of Formation of Weakly Ferromagnetic Strings (PFWFS) and that we are going to 
calculate this probability in the saddle point approximation. Therefore, we consider a configuration 
of the field (instanton) which satisfies the boundary condition imposed by the PFWFS (|91l92ft 

dt#(x,t)\ t=0t0<x<n = p (CI) 

and that minimizes the action, i.e. that satisfies the Euler-Lagrange equations: 

(dpd? - m 2 ) = 0. (C2) 

The latter equation is the Klein-Gordon equation with the mass given by m 2 = 2 7 (see ©). The 
PFWFS will be found from the value of the action So corresponding to this instanton configuration 
(|94|) , In this appendix we calculate the stationary action needed in Sec. IVIlTI 

We now solve the differential equation (|C2|) with non-trivial boundary condition (|C1|) by re- 
casting it as the integral equation: 

1 rn . . 

0OM) = ^J dtK ° (WO* - v) 2 + 1 2 ) f(y) dy, (C3) 

where Ko(x,x';t,t') is the modified Bessel function of 0-th order - the kernel of the differential 
operator ()C2|) in two dimensions. We impose the boundary condition (|C1|) by requiring that the 
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1 r { „ / n ^ ^ m 2 t 2 



"source" f(y) satisfies 

dt$(x,0)\ Q<x<n = limi- [jif 2 (m^-tf + i 2 ) 



(x - y y + 1 2 

(mV(x-y)* + t*) V(x _ m y)2 + t2 ) M *y = p- (C4) 



This is the integral equation on f(y) we have to solve. 

Once the limit t — > is taken, the kernel in Eq. (|C4j) is singular. We isolate the singularity by 
rewriting equation (|C4|) as: 

d r m_ dy + lim r Go(X5 1 )/( j dy = 27r - (C5) 

dx y x - y t->Q J 

with 

(■£ — — ^2 , . Tn 2 t 2 

G «(^y1 s ( X _ y) 2 + t 2 + g 2 ( W(* - yf + t'j (i _ y? + <2 

-^( m ^^)7(=FTP (C6) 

or, after integration over x, as 

T iM. d y + r G(x; d y = 2 TT f) X (C7) 

Jo x-y Jo 

with 

G{x-y) = \im! Go(x x ,t;y) da*. (C8) 
Jo 

We have recasted Eq. ()C4|) in the standard form for a singular integral equation l)C7|) . Once we 
have the solution of this equation, we can calculate the action corresponding to this instanton as 



So = 4V1 - /i 2 P / /(y)dy. (C9) 



o 

We solved singular integral equation (|C7|) numerically and we computed the corresponding 
action l)C9|) . The results of these calculation are presented as a plot of the action S$ vs. n in 
Fig. [3 where we notice the crossover from a quadratic to a linear behavior (corresponding to a 
crossover from Gaussian to exponential behavior for the probability, (|94[)) as we expected. To 
confirm the nature of this crossover, in Fig. Elwe plot dSo/dn and we see that it starts linearly and 
then saturates asymptotically as it should. 

In the limit n <C 1/m, we can expand the Bessel functions in the kernel QC6|) 

m 2 ( t 2 1 , r , x9 01 , m „ 1 



Go(x, t; y) = -— ( (a; _ y) 2 + t 2 + 2 In [(x - y) 2 + t 2 ] + In y + G - - ) + . . . , (CIO) 
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FIG. 9: The solid line is the plot of the stationary action l|C14|) against n. This analytical solution is valid 
for n « 1/m and corresponds to m = 0.01 and p = 0.2. The dotted line represents the value of the 
action (|G9|) with the source given by numerical solution of the singular integral equation (|C7|I . The dashed 
line corresponds to the zeroth-order, pure Gaussian, solution, i.e. 1C14I) with m = 0, which we include 
for comparison. We see that the inclusion of first order correction almost doubles the range in which the 
analytical solution is accurate. 



where G is Catalan's constant. Then we solve the singular integral equation QC7|) to first order by 
first transforming it into a regular integral equation. 

In 3l|, Chap. 14, Sec. 114 it is explained that a singular integral equation like (j(JT|) is equivalent 



to 



f(x) + — / N(x; y) f(y) dy = 2py/x{n - x), 
m Jo 



(Cll) 



where the new kernel is 



N(x;y) 



^/x{n - x) f n G(y';y) 



£ 



dy'. 



(C12) 



vri Jo y/ y '(n-y')(x-y') 

Using (|C10|) . we can explicitly calculate the integral defining N(x;y) in terms of elementary 
functions and after some algebra the integral equation (|C11|) results in a long, but essentially 
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simple, regular integral equation. Its solution is 

f(x) = p 



2 + 



2 2 

m n 



, mn _ 3 

In hG 

8 2 



\/ x(n — x) — p 



2 2 

m n 



n 

x tan 

2 



-i 



n — x 



The corresponding stationary action (|C9j) is 



So = 7rV 1 - h 2 p 



2 2 

n 



1 + 



2 2 

m n 



16 



i run „ _ 
In hG-2 



(C13) 



(C14) 



The first term in ()C14|) corresponds to the Gaussian decay of PFWFS we expect in the limit of 
m = 0. In Fig. |H1 we compare this analytical result for the action with the numerical result of 
FigEl In the plot, we include the pure Gaussian decay (the first term in ()C14jO . which already 
gives a remarkable agreement for small n. The full solution ()G14|) extends this agreement further 
for larger n. 
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